Draft version January 12, 2013 

Preprint typeset using I^'T^]X style emulatcapj v. 11/10/09 



O 



A DEEP VIEW OF THE MONOCEROS RING IN THE ANTICENTER DIRECTION: CLUES OF ITS 

EXTRA-GALACTIC ORIGIN 

A. SoLLiMA^'^'*, D. Valls-Gabaud'', D. Martinez-Delgado'"'^'^, J. Fliri^, J. Penarrubia"^' , H. Hoekstra'^'** 

Draft version January 12, 2013 

ABSTRACT 

We present the results of deep imaging obtained at the CFHT with MegaCam in the Anticenter 
direction at two different heights above the Galactic disk. We detect the presence of the Monoceros ring 
in both fields as a conspicuous and narrow Main Sequence feature which dominates star counts over 
a large portion of the color-magnitude diagram down to g' ^ 24. The comparison of the morphology 
and density of this feature with a large variety of Galactic models excludes the possibility that it can 
be due to a flare of the Galactic disk, supporting an extra-Galactic origin for this ring-like structure. 
Subject headings: Galaxy: structure — galaxies: interactions — Galaxy: disk — methods: data 
analysis — techniques: photometric 



< 

o 

o 



> 

O 



1. INTRODUCTION 

The Monoceros ring represents one of the most contro- 
versial features among the various sub-structures discov- 
ered in the Milky Way. It was discovered as a ring-like 
structure spanning about 170° in longitude at a nearly 
constant Galactocentric distance at low-latitudes (New- 
berg et al. 2002; Yanny et al. 2003). Its structure and 
extent have been later studied by many authors using 
both photometric (Rocha-Pinto et al. 2003; Ibata et 
al. 2003; Conn et al. 2005a, 2007) and spectroscopic 
(Crane et al. 2003) data. On the basis of these anal- 
yses, the Monoceros (Mon) stream has been found on 
both sides of the plane of the Galaxy at Galactocentric 
distances 15 < i? < 20 kpc. Subsequent spectroscopic 
studies indicated that this structure is kinematically cold 
(13 < cr„ < 20 Km s-i; Conn et al. 2005b; Martin et 
al. 2006; Casetti-Dinescu et al. 2008) and homogeneous 
in metal content {<J\FelH\ ~ 0-15 dGX\ Ivezic et al. 2008, 
hereafter 108). The nature of the Mon ring is highly de- 
bated: in fact it is not clear if this feature is a debris from 
an accreted satellite (e.g. Helmi et al. 2003; Martin et 
al. 2004; Penarrubia et al. 2005; Martinez-Delgado et al. 
2005) or an intrinsic feature of the Galactic disk possibly 
associated to the disk flare (Momany et al. 2006, here- 
after M06; Hammersley & Lopez- Corredoira 2011), the 
Norma-Cygnus spiral arm (Moitinho et al. 2006), pertur- 
bations due to past accretion events (Kazantzidis et al. 
2008) or the second dark matter caustic ring (Natarajan 
& Sikivie 2007). 
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In this Letter we present the detection of the Mon ring 
toward the Galactic anticenter at two different Galactic 
latitudes and use its color-magnitude diagrams (CMDs) 
to investigate the possible connection between this ob- 
ject and the Galactic disk flare. The Anticenter rep- 
resents a particularly convenient direction to test this 
scenario since it lies close to the line of nodes where 
the effects of the Galactic warp are minimized (Drimmel 
& Spergel 2001), ii) it avoids contamination from bulge 
stars, and Hi) it lies along a radial away from the Galac- 
tic center, following the direction of growth of the disk 
scale-height. This analysis comes from two photometric 
campaigns performed with MegaCam@CFHT devoted to 
the study of weak lensing in massive clusters of galax- 
ies (see Hoekstra 2007) and to the search for extra-tidal 
structures in the outskirts of outer halo globular clusters 
(Martinez-Delgado et al. 2004). 

2. OBSERVATIONS AND DATA REDUCTION 

The photometric data have been obtained with the 
MegaCam camera, mounted at the Canada-France- 
Hawaai Telescope (CFHT)(Mauna Kea, Hawaai). The 
camera consists of a mosaic of 36 chips with a pixel 
scale of 0.185" pixel~^ providing a global field of view 
of ^ 1° X 1°. Observations have been performed in 
Service Mode in two different observing runs on Octo- 
ber 2005 (3 nights) and January 2009 (two nights). The 
two datasets consists in a set of 15 (4 g' and 11 r') 600 
sec-long and 13 (6 g' and 7 r') 680 sec-long images cen- 
tered at {a,S) = {7h 18m 08s, +37° 37' 45")(field 1) 
and {a, 5) = (Ih 38m 08s, -1-38° 54' 00")(field 2) with a 
dithering pattern of few arcminutes to fill the gaps be- 
tween the chips. The average seeing for the two datasets 
was 0.7" and 1.0", respectively. 

The standard reduction steps (bias, dark and flat-field 
correction) have been performed using the Elixir pipeline 
developed by the CFHT team. We used DAOPHOT II 
and the PSF fitting algorithm ALLSTAR (Stetson, 1987) 
to obtain instrumental magnitudes for all the stars de- 
tected in each frame. Mean frames have been obtained 
by aligning and averaging all the g' and r' images with a 3 
a clipping rejection thresold, using the TerapisQ software 

^ http://tcrapix.iap.fr 
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specifically devoted to this task. The automatic detec- 
tion of sources has been performed on the mean frame 
adopting a 3a treshold. The file with the object positions 
has been then used as input for the PSF-fitting proce- 
dure, that has been performed independently on each im- 
age. As usual, the most isolated and bright stars in each 
field have been employed to build the PSF model (here 
a Moffat function of exponent 2.5 has been adopted). 
For each passband, the derived magnitudes have been 
transformed into the same instrumental scale and aver- 
aged. We adopted the nightly zero points and reddening 
coefficients provided by the CFHT team to link the in- 
strumental magnitudes to the standard system. Finally, 
a catalog with more than 40,000 calibrated sources has 
been produced. The photometrically calibrated catalog 
has been cross-correlated with the seventh data release 
of the Sloan Digital Sky Survey (SDSS; Abazajian et 
al. 2009) catalog, which lists accurate positions for more 
than 3,000 objects over an area of tt sq.deg. around both 
fields, to derive an accurate astrometric solution with a 
typical r.m.s of 200 mas. These catalogs represent the 
deepest and most accurate datasets on the Mon feature, 
with over 3 times more stars and a photometric accuracy 
about an order of magnitude better than previous SDSS 
data. 

3. RESULTS 
3.1. Color-Magnitude Diagrams 

In Fig. [U the g' - r' , g' CMDs of the two fields 
are shown. Only objects with a sharpness parameter 
\S\ <0.5 (as defined by Stetson 1987) have been plotted 
to avoid contamination from background galaxies. As 
can be noted, the CMDs sample the foreground Galac- 
tic population, reaching a limiting magnitude g' 25.5. 
A significant overdensity of stars can be noticed in both 
fields at g' — r' < 1 and 19 < 5' < 24, with a morphol- 
ogy which resembles a Main Sequence (MS). It is evident 
from Fig. [1] that the density of MS stars is significantly 
smaller in field 2 compared with field 1. 

The present detections are located close to the re- 
gion where the Mon ring was discovered (Newberg et al. 
2002). More recently, Grillmair (2006) found a compact 
overdensity of stars above the background field popula- 
tion at an height of ~30 degrees above the Galactic plane 
close to the previous detection of Mon (the "Anticenter 
stream"). Although the possible association between the 
Anticenter stream and Mon is still debated (see Carlin et 
al. 2010 and references therein), our target fields are lo- 
cated at lower Galactic latitudes (A5 ~ 12°) with respect 
to the Anticenter stream overdensity, being therefore un- 
related with such structure. 

3.2. Synthetic CMDs 

To test the possible association of the MS feature ob- 
served in Fig. [T]with the disk flare we simulated a series 
of synthetic CMDs for the Galactic field population fol- 
lowing the prescription of seven different Galactic mod- 
els: Drimmel & Spergel (2001; hereafter DSOl); Lopez- 
Corredoira et al. (2002; L02); Robin et al. (2003; R03), 
Foster & Routledge (2003; FR03), Yusifov (2004; Y04), 
Hammersley & Lopez-Corredoira (2010; HLll) and a 
simple model which has been calibrated to fit the 2MASS 
data of the disk flare and warp collected by M06. Here- 
after, we will refer to this best-fit model as Sll. For this 
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Fig. 1.— CMDs of field 1 (left panel) and field 2 (right panel). 
Only stars with a sharpness parameter \S\ < 0.5 are plotted. The 
clump visible at g' > 24, g' — r' <0.5 is due to the contamination 
of unresolved background galaxies. 



TABLE 1 

Best-fit parameters of the best-fit Sll model 



Parameter Value Description 

Rq 8.5 kpc Solar Galactocentric distance 

Zq 15 pc Solar height above the Galactic plane 

g 0.63 Halo flatness parameter 

/ 0.9 Fraction of thin disk stars 

0.001 Fraction of halo stars 

2.77 Halo density power-law coefficient 

R\ 1.99 kpc Thin disk scale-length 

R.2 2.99 kpc Thick disk scale-length 

hz,i(RQ) 193 pc Thin disk scale- height at the Solar circle 

'^2.2(^0) 611 pc Thick disk scale-height at the Solar circle 

hji 23.0 kpc Flare characteristic radius 

Zo 370 pc Maximum warp height 

Rw 11.5 kpc Warp starting radius 

3.0 kpc Warp characteristic radius 

00 11° Line of nodes angle 



last model we adopted an oblate spheroidal halo with a 
declining profile of the form 

Ph{R, Z) = h PQ,d (rq/VWTWWY' (1) 

The disk has been modelled as a sum of two exponential 
disks (Juric et al. 2008; hereafter J08) 

Pd{R, Z) - po^, (/ e-(«-«o)/«i-(^+^e)/^-.i + (2) 

(1 _ y) g-(i?.--Ro)/i?.2-(2+2o)//iz,2) 

The scale-heights of both disks have been assumed to 
increase with distance from the Galactic center to take 
into account the disk flare 

/iz,z(i?)=/i2.^(i?0) e(^-^«)/''^' (3) 



The Monoceros ring in the Anticenter 
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Fig. 2. — Disk scale- height as a function of the Galactocentric 
distance. The measures by M06 are plotted as filled circles. The 
prediction of the various models (DSOl: blue dot-dashed line; L02: 
red solid line; R03: red dot-dashed line; FR03: blue solid line; 
Y04: blue dashed line; Sll: black solid line; HLll: red dashed 
line; "wall flare" model; black dashed line) arc overplotted. 

Finally, the disk warp has been modelled as a series of 
tilted rings whose disk heigths is defined as 

Zj^{Z,(f>) = Zq tanh{{R— Rw) / ht^,) sm(0— ^o) for R > R., 

.(4) 

where (f> is the angle with respect to the Sun direction 
in a cylindrical system of coordinates centered on the 
Galactic center. In the definition of our best-fit model 
we kept fixed those parameters derived by JOS and left 
free the 4 ones which determine the warp and flare shape 
{hfi, Zq, hw and (po). The final values of all the param- 
eters involved in the above model are listed in Table 1. 
The value of pQ^d has been normalized to reproduce the 
observed number of stars in field 1. 

All the above models take into account for both the 
disk flare and warp and have been calibrated to repro- 
duce the density of different markers: warm dust (DSOl), 
bright stars (L02), dwarf stars (R03, HLll), HI gas 
(FR03), pulsars (Y04) and red clump stars (M06). In 
Fig. [2] the prediction of the various models considered in 
this work arc compared with the data collected by M06. 
Note that signiflcant differences are evident between the 
various models mainly resulting from the different ob- 
jects used as disk markers. 

The overall CMD has been simulated as the sum of 
the contribution of the halo and disk stellar populations. 
We ignored the bulge contribution assuming it negligible 
in the Anticenter direction. For each Galactic compo- 
nent we randomly extracted a population of stars from 
a Kroupa (2001) mass function with suitable ages and 
metallicitiej*^. Only R03 and HLll distinguish between 

Being the simulated CMDs mainly constituted by unevolved 
stars, they are not affected by the uncertainties on the adopted SFR 
for the various Galactic components, but mainly on their metal- 



thin and thick disk adopting two different characteristic 
lengths and heights. For the other models we simulated 
the Galactic disk population as a mixture of thin and 
thick disk stars in proportions estabilished by JOS. This 
approximation ensures a good level of accuracy in re- 
producing the CMD (see also DSOl; L02; FRQ3; MQ6; 
Reyle et al. 2009). We adopted a two-step star for- 
mation rate (SFR) and the age-metallicity relation from 
Fuhrmann (199S) for the thin disk (Girardi et al. 2005), 
while for the thick disk a constant SFR between 11 
and 12 Gyr and a gaussian metallicity distribution with 
[M/H] = —0.4 ± 0.1 has been assumed (Girardi et al. 
2005). For the halo, we assumed a constant SFR be- 
tween 12 and 13 Gyr and the metallicity distribution de- 
rived by Ryan & Norris (1991). A fraction of binaries of 
/ = 50% (Duquennoy & Mayor 1991) has been simulated 
by randomly associating pairs of stars and calculating the 
corresponding magnitudes by summing the fluxes of the 
two components (see Rubenstein & Bailyn 1997). Abso- 
lute magnitudes and dereddened colors have been then 
derived using the evolutionary tracks by Marigo et al. 
(2008) in the MegaCam photometric system. Halo and 
disk stars have been then located at different distances 
according to the prescription of each Galaxy model. For 
each star a proper extinction has been included by cal- 
culating the dust column density at the star's distance 
assuming a dust distribution across the disk of the form 
Pdust tx e~^^^'^"'*sech?{Z/Zdust) with a scale- heigth of 
Zdust — 134.4 pc and a scale-length of Rdust =2.26 kpc 
(DSOl), and normalized by imposing the extinction at in- 
ilnity predicted by the Schlegel et al. (1998) maps. The 
distance and extinction have been finally used to convert 
absolute magnitudes and colors into apparent ones. Pho- 
tometric errors and incompleteness have been estimated 
by means of artificial stars experiments which have been 
performed on the CFHT images. A detailed description 
of the adopted procedure can be found in Bellazzini et 
al. (2002). 

In Fig [3] the observed CMD of field 1 and the cor- 
responding synthetic CMDs are shown. It is evident 
that while models DSOl, R03, FR03 and Sll wefl re- 
produce the overall morphology of the Galactic field evo- 
lutionary sequences, none of them is able to reproduce 
the prominent MS feature clearly visible in the observed 
CMD. A poor representation of both the overall CMD 
and the Mon MS is instead given by Y04 and HLll mod- 
els. Only model L02 predicts a compact distribution of 
stars resembling a MS. Such an increase is due to the 
extreme fiaring of the disk adopted by L02 (see Fig. [2]) 
which produces a sudden increase of the density at a 
well defined Galactocentric distance (R~ 14 kpc), while 
the double-exponential density law [pd oc ^-R/hR-z/hz^^ 
dampen the distribution of stars at longer distances. 
Note however that this model (which has been extrapo- 
lated here outside its range of validity; i? < 15 kpc) fails 
to reproduce the magnitude spread of the observed MS 
feature as well as the morphology of both the red (at 
g' < 18.5, g' -r' < 0.7) and blue (at g' - r' > 1) field 
dwarfs in Fig. [T] 

Summarizing, it seems impossible to reproduce the ob- 
served morphology of the CMDs with the existing models 

licity distributions, which are well constrained by spectroscopic 
observations (Venn et al. 2004). 
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Fig. 3. — Comparison between the CMDs of field 1 and the seven Galactic models considered in this work. 



of disk flare. To check the general validity of such con- 
clusion we compared the observed CMD with a model 
where the most extreme disk flaring has been assumed. 
For this purpose we replaced in model SIX the functional 
form of the disks flaring (eq. [S]) with a "wall" function 
where the scale-height of the two disks abruptly goes to 
infinity at a given Galactocentric distance (set to 15.5 
kpc to fit the mean magnitude of the Mon MS; see Fig. 
[2]). The simulated CMD for this model (hereafter re- 
ferred as "wall flare") is shown in the central panel of 
Fig. m In Fig. |4]we compared the distribution of mag- 
nitude differences about the MS ridge line in the color 
interval 0.4 < g' — r' < 0.9. It is apparent that the 
"wall flare" model predicts a width of the MS which is 
not compatible with the observed one: the standard de- 
viation of the distribution is (J^g' — 0.30 ± 0.03 in the 
observed CMD and aAg' = 0.85 ±0.08 in the "wall flare" 
model (i.e. ^ 3 times larger). Only an "ad hoc" abrupt 
cutoff of the density profile of both disks at a distance 
close to the flare radius (in contrast with what found by 



J08) would reduce the magnitude spread of this model. 
A second test has been performed by comparing the ratio 
between the number of MS stars observed in field 1 and 2 
with the model prediction. For this purpose, we defined 
a region in the CMD containing all stars with g' mag- 
nitude within 0.5 mag about the MS ridge line defined 
in field 1, a color 0.4 < g' — r' < 0.9 and a magnitude 
g' < 23.8, to avoid contamination from Galatic dwarf 
stars and to ensure a level of completeness V' > 90% 
in both fields. While the fraction of MS stars appears 
to increase by a factor N1/N2 — 1.87±0.11, the "wall 
flare" model predicts a small decrease of this fraction 
(A^i/iVa = 0.92 ± 0.03 JB- Therefore, on the basis of the 
above analysis we conclude that even the most extremely 

The N1/N2 population ratio calculated in the "wall flare" 
model turns out to be smaller than unity since the small shift in 
magnitude produced by the decrease of latitude moves part of the 
MS feature outside the adopted selction box. This is a consequence 
of the large magnitude spread of the MS feature predicted by this 
model. 




Fig. 4. — Upper panels: Zoomed CMDs of field 1 (upper-left panel), "wall flare" model (upper-central panel) and simulated stream model 
(upper-right panel). The MS ridge line and the color range used to construct the magnitude distribution are shown with red solid and 
dashed lines, respectively. Bottom panel: Distribution of magnitude differences about the MS ridge line for the observed CMD of field 1 
(black histograms), "wall flare" model (red histograms) and the simulation for a discrete stream (blue histograms). The best-fit gaussians 
are overplotted. 



acceptable flared disk model cannot reproduce the ob- 
served CMD since it fails to predict: i) the population 
of both blue and red field dwarfs; ii) the 2MASS data of 
the flare measured by M06 (see Fig 12]), Hi) the magnitude 
dispersion of the observed MS feature, and iv) the ratio 
of MS stars at the two observed Galactic latitudes. 

For comparison with the stream scenario, we con- 
structed the synthetic CMD of a galaxy remnant located 
at a given distance to the Sun with a variable thickness. 
For this purpose, we apply the Bayesian formalism de- 
veloped by Hernandez & Valls-Gabaud (2008) to infer 
the posterior probability distribution functions of dis- 
tance and age, using as priors a gaussian function for 
metallicity, centred on [Fe/H] — —0.95 and with disper- 
sion 0.15 (108). We obtain maxima in the marginalised 
posterior probability distribution functions at an age of 
t = 9.2 ± 0.2 Gyr and a distance oi d = 9.1 ± 0.2 kpc 



with a dispersion of Od — 0.90 ±0.08 kpc, consistent with 
previous estimates for the Mon stream (Newberg et al. 
2002) and the predictions of Penarrubia et al. (2005). 
The obtained CMD is shown in the right panel of Fig. 
01 In this last case the agreement with the observed 
CMD is striking. A good agreement can be also found 
by assuming the Mon ring formed by stars with the same 
age and metallicity distribution of the Galactic disk (e.g. 
the case of a disk perturbation or a spiral arm). In this 
case, the best-fit distance and dispersion turn out to be 
d — 12.6 ± 0.8 kpc and ad — 0.3 ± 0.2 kpc, respectively. 
However, it is worth noting that while agreement can be 
had over the limited magnitude range used for estimating 
distance, the metallicity estimates by 108 argue against 
this scenario. 

4. DISCUSSION 
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In this Letter we present the results of a deep Mega- 
Cam imaging of two regions in the Galactic Anticenter 
direction where the presence of the Mon ring is clearly 
evident. The high quality of our CMDs allowed to com- 
pare the morphology of the observed MS feature with the 
predictions of the most recent Galactic models including 
the effect of disk flare and warp. We find that none of the 
considered models is able to reproduce the position and 
morphology of the Mon MS or the ratio between the num- 
ber of MS stars in the two observed fields, even with ex- 
treme fiares. We therefore conclude that, in the absence 
of an " ad hoc" abrupt disk cutoff, the photometric signa- 
ture of the Mon ring cannot be explained by any smooth 
variation of the Galactic disk structure. This result is in 
contrast with the conclusion drawn by HLll. The rea- 
son of such a discrepancy stems from the method used 
by these authors to calculate the predicted star counts 
of their model: indeed, they adopted a delta function in 
magnitude to calculate the number of stars displaced at 
a given distance to the observer neglecting many impor- 
tant factors like the slope of MS stars in the CMD, the 
composite stellar population of the Galactic components 
and the photometric errors that all contribute to increase 
the magnitude spread of MS stars. By neglecting these 
effects HLll interpreted the overall magnitude spread of 
the Mon MS as entirely due to a distance spread, over- 
estimating its actual dispersion. 

Of course, in absence of direct spectroscopic data, our 
data cannot exclude the other possible interpretations 
of the Mon ring as the result of small scale perturba- 
tions of the Galactic disk due to various possible physi- 
cal processes (see e.g. Moitinho et al. 2006; Natarajan 
& Sikivie 2007; Kazantzidis et al. 2008). However, our 
CMDs impose strong constraints on the vertical structure 
and extent of such a structure. Indeed, if we assume the 
Mon ring formed by Galactic disk stars, the location in 



the CMD of the observed MS feature and its dispersion 
in magnitude are compatible with an object located at 
a distance d = 12.6 ± 0.8 kpc to the Sun and a width 
of FWHM < 1 kpc. This corresponds to a very com- 
pact structure located at an height above the Galactic 
plane oi Z ^ 4.5 kpc (i.e. ~ 4.5 times the local scale- 
height of the disk). Moreover, the density of Mon stars 
decreases of almost a factor of two by moving only 4° 
toward higher Galactic latitudes (see Sect. 13.21) . On the 
other hand, under the assumption that the Mon ring is 
due to an accretion event which occurred in a past epoch, 
the derived distance of c? = 9.1 ± 0.2 kpc and dispersion 
fJd = 0.90 ±0.08 kpc (FWHM=2.1±0.2 kpc) are in good 
agreement with the prediction of the model by Penarru- 
bia et al. (2005; see their Fig. 7). 

Our results lend support to the previous studies on 
the metallicity distribution (108), abundance patterns 
anomalies (Chou et al. 2010) and kinematics (Casetti- 
Dinescu et al. 2008) of this object which show that the 
stellar composition of the Mon ring is different from that 
of the Milky Way, favouring the extra- Galactic scenario. 
An accurate high- resolution spectroscopic survey of these 
faint MS stars would defitively clarify the real nature of 
the Mon ring. 

Based on observations obtained with MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the CFHT ob- 
serving programs 05BC19 and 09AF03. This work is 
based on data products produced at TERAPIX and the 
Canadian Astronomy Data Centre. This research was 
supported by the Spanish Ministry of Science and In- 
novation (AYA2007-65090), by ANR POMMME (ANR 
09-BLAN-0228) and CNRS/MAE PICASSO. We warmly 
thank the anonymous referee and M. Lopez-Corredoira 
for their helpful suggestions and Y. Momany for provid- 
ing his data. 
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